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ABSTRACT 



We have performed RV monitoring of the components of the binary system HD 106515 over about 11 years using the 
high resolution spectrograph SARG at TNG. The primary shows long-period radial velocity variations that indicate 
the presence of a low mass companion whose projected mass is in the planetary regime (msini = 9.33 Mj). The 9.8 
years orbit results quite eccentric (e = 0.57), as typical for massive giant planets. Our results confirm the preliminary 
announcement of the planet included in Mayor et al. (2011). The secondary instead does not show significant RV 
variations. The two components do not differ significantly in chemical composition, as found for other pairs for which 
one component hosts giant planets. Adaptive optics images obtained with AdOptQTNG do not reveal additional stellar 
companions. From the analysis of the relative astrometry of the components of the wide pair we put an upper limit 
on the mass of the newly detected companion of about 0.25 Mq. State of art or near future instrumentation can 
provide true mass determination, thanks to the availability of the wide companion HD106515B as reference. Therefore, 
HD106515Ab will allow deeper insight in the transition region between planets and brown dwarfs. 
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1. Introduction 

The upper mass limit of planetary objects is currently 
widely debated in the scientific community. On one hand, 
an operational definition can be based on the minimum 
mass for deuterium burning (about 13 M j) as a divid- 
ing li ne between planets and brown d warfs ([Burrows et al.l 
120011) . as defined by IAU (|Boss et al.l feoij. orJo_differ- 
ent fixed threshold values (e.g. 24 — 25 Mj: iButler et al.l 
120061: iSchneider et~ai1l2011l) . On the other hand, a defini- 
tion based on formation mechanisms is more difficult to 
obtain. In fact, observational data are incomplete and the 
detection techniques in most cases can provide only indi- 
rect indication on the formation of a detected substellar 
objects. Furthermore, only minimum masses are known for 
most of the objects detected by RV surveys or, in the case 
of objects detected through direct imaging, true masses are 
uncertain because of the large sensitivity of luminosity on 
age and intrinsic un certainties of theore tical models espe- 
cially at young ages (jBaraffe et al.ll2003t) . 



Send offprint requests to: S. Desidera, 
e-mail: silvano . desideraOoapd. inaf . it 

* Based on observations made with the Italian Telescopio 
Nazionale Galileo (TNG) operated on the island of La Palma by 
the Fundacion Galileo Galilei of the INAF (Istituto Nazionale 
di Astrofisica) at the Spanish Observatorio del Roque de los 
Muchachos of the Instituto de Astrofisica de Canarias. 



In spite of these difficulties there is growing evidence 
for a significant overlap in mass between objects formed 
like stars do and objects formed in a protoplanetary disks. 
Several objects of planetary mass have been detected as 
free floating objects in star cluster s, star forming regions 
or in the field using imaging (e.g. IZapatero Osorio et al.l 
l2000l:lLeggett et al.l[20lllScholz et alJ2012f ) and microlens- 
ing (ISumi et al.l 1201 lh or as very wide companions of stars 
(e.g. IChauvin et aLll2005l) . Some of these objects might be 
formed in planetary systems and then pushed at very wide 
separation or ejected from the system because of dynami- 
cal in teractions with other planets (Ju mping Jupiters sce- 
nario, (MarzarT&|Weidcnschilling 2002). However, it seems 
unlikely this is the only mechanism produci ng free float- 
ing o bjects below deuterium burning mass (jBowler et al.l 
2011). Instead, the minimum mass for core collapse was 
found to be of a few Mj and is then likely that objects 
of planetary mass form ed star-like outside planetary disks 
(jWhitworth et al.ll2007ri . The statistics of low mass brown 
dwarfs appear to be compatible with Je ans mass fragmenta- 
tion of an interstellar molecular cloud (jZuckerman fc Sond 
l2009h . 

On the other hand, there are objects with masses from 
13 to about 25 Mj that are found in systems with other 
lower mass planets, such as HD168443 and HAT-P13. In the 
latter case the planetary nature of the lower m ass compan- 
ion is confirmed by the occurrence of transits (jBakos et al.l 
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2009). In other cases such as HD 38519 and HD 202206 a 
debris disk is present in t he system beside a massiv e planet 
and a lower mass planet ()Moro-Martrn et aLll2010f) . These 
facts support the formation of these objects in a proto- 
planetary disks. Planet formation models also predict the 
presence of very massive planets, up to about 38 Mj, in ex- 
ceptional cases of long- lived, massive and metal-rich disks 
(|Mordasini et al.l 120091). The rising mass function below 



about 20 — 30 Mj (jGrether fe Lineweaverll2006f ) is another 
indication that a different formation mechanism starts to 
be present above deuterium burning mass. 

This picture is further complicated by evidences that 
the statistical properties of planetary mass companions 
with (projected) masses between 4 to 15 Mj are different 
from those of lower mass planets (|Ribas fe Miralda-Escudel 
l2007h . It is then possible that either a different formation 
mechanism is in action or that the evolution of massive 
planets in a protoplanetary disk is different depending on 
planetary mass. 

The statistic of substellar objects in the mass range be- 
tween 10 to 30 Mj i s still limited, given the intrinsic rar- 
ity o f these objects (jSozzetti fc Desidera] l2010t iDiaz et all 
|2012| ). The discovery of additional candidates is then wel- 
come, especially when the true mass of the companion can 
be de termined or significant l y constrained through astrom - 
etry (jSahlmann et all 120111: iReffert fc Quirrenbacb]|2011l) . 
We present here the confirmation of a high-mass planet 
candidate orbiting the star HD106515A and clues on its 
mass from astr ometry. The planet was first included in the 
compilation bv lMavor et al.l (|201lf) but only orbital period, 
RV semiamplitude, eccentricity, and corresponding mini- 
mum masses and semimajor axis are listed, postponing a 
more detailed analysis to a forthcoming paper. HD 106515A 
is part of a wide binary system, with the companion HD 
106515B at a projected separation of about 250 AU. Both 
components were observed as part of the RV survey look- 
ing for planets around the components of moderately wide 
binaries performed using SARG at TNG. 



2. Observations and data reduction 

Observations were performed at the Italian Telescopio 
Nazionale Galiel eo (TNG) using the high resolution spec- 
trograph SARG (jGratton et al.lhOQll ). All but one spectra, 
used as template in the RV determination process, were 
acquired with the iodine cell inserted in the optical path. 
The observing procedure and the instrument set-up are the 
st andard ones for SAR G planet search survey. We refer 
to iDesidera et al.l (|201ll ) for further details. In the case of 
HD106515, the chance of contamination of the spectra is 
negligible even in the worst observing conditions, thanks to 
the 6.8 arcsec separation on the sky between the compo- 
nents. 

The acquisition of data on HD106515 was stopped after 
May 2009 because of lack of time all ocation and was re - 
cently restarted after the publication of lMavor et al.l (|2011l) , 
that listed a planet around the primary component of this 
system. The 900s integration yields typical S/N ratios of 
100. for HD106515A and 90 for HD106515B. 

Differential ra dial velocities w ere obtained using the 
AUSTRAL code (jEndl et al.ll2000D . achieving internal er- 
rors of about 3 m/s. Data taken in 2011-2012 have a some- 
what larger uncertainty, because of significant asymmetries 



of the spectrograph instrument profile. Standard RV anal- 
ysis performed using a central Gaussian with two or four 
Gaussian satellites show RV higher by few tens of m/s with 
respect to previous data. However, when analyzing the data 
exploiting the Maximum Entropy algorithm such discrep- 
ancy vanishes. In any case, while the internal errors are 
similar to the older data, we can not exclude the occurrence 
of systematics at 10-15 m/s level for these recent data. The 
RV time series of HD106515A and B are reported in Tabled 
and[U 

HD 106515 was also observed o n 21 June 2007 using 
AdOpt@TNG (|Cecconi et al.ll2006l) . At that time only a 
long term trend was appearing from RV data and our ob- 
servations were aimed at the direct detection of the com- 
panion responsible for the trend. The acquisition and data 
re duction procedures ar e the same as done for HD132563 
in IDesidera et al.l (|2011h , with the differences that all the 
images on HD106515 were taken using the broad band K' 
filter and at the same rotation angle. 194 images were ac- 
quired and used in our analysis. The projected separation 
between HD106515 A and B results of 6.897±0.015 arcsec 
and the position angle of 267.07±0.12 deg. Detection limits 
in K band magnitude difference were transf ormed in mass 
limits using the mass-luminosity relation by iDelfosse et al.l 
(I2000D at Mk < 9.5 a nd the 5 Gyr theoretical isochrone by 
IChabrier et al.l (|2000D at fainter magnitudes. 



3. Stellar properties 

HD 106515 (HIP 59743, GJ 9398, ADS 8477) is a pair 
formed by two similar stars slightly less massive than the 
Sun. The main properties of the components of HD 106515 
are summarized in Table [T] 

The metallicity of the components was studied by us us- 
ing bo th SARG and FEROS spectra (jDesidera et al.ll2004 
l2006al ). The metallicity resulted close to solar. The line- by- 
line abundance differential analysis did not reveal signifi- 
cant abundance difference between the components. 

The slow projec ted rotational velocity , low level of chro- 
mospheric activ ity ((Desidera et al.|[2006at iGrav et al.ll2003l : 
lArriagadall201lD. and lack of detection in ROSAT All Sky 
Survey (|Voges "etan[2000l) are all consistent with a rather 
old age, at least as old as the Sun. Isochrone fitting does 
not allow to put firm constraints on stellar age, because of 
long timescales of the evolution of stars less massive than 
the Sun. The thin disk kinematic places an upper limit of 
about 8 Gyr. The system age is then in the range between 
4 to 8 Gyr. 



4. A massive planet around HD106515 

The RV time-serie of HD 106515A shows a long period mod- 
ulation clearly detectable by eye (Fig.rjJ). Lomb-Scargle pe- 
riodogram confirms the long periodicity and the lack of 
shorter periodicities in the data. The low level of chro- 
mospheric activity guarantees the Keplerian origin of the 
RV variations observed for HD 106515A. A n expec t ed ac- 
tivity jitter of 3.5 m/s was derived using Wright] (|2005h 
calibration. We then performed Keplerian fitting using 
a Levenberg-Marqua rdt least-squares fit algorithm as in 
IDesidera et al.l (1201 if). The pa rameters are similar to those 
derived bv lMavor et al.l (|2011l ). Both orbital parameters are 
listed in Table [3J The residuals from the best fit Keplerian 
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Table 1. Stellar properties of the components of HD 
106515. 



Parameter HD 106515 A HD 106515 B Ref. 



a (2000) 12 15 06.567 12 15 06.103 1 

5 (2000) -07 15 26.38 -07 15 26.61 1 

fi a (mas/yr) -249.67 ± 0.91 2 

Us (mas/yr) -52.29 ± 0.74 2 

RV (km/s) 20.66± 0.11 19.94±0.11 3 

7T (mas) 28.42 ± 0.96 2 

d (pc) 35.2±1.1 2 

U (km/s) -28.0 4 

V (km/s) -38.7 4 
W (km/s) 4.6 4 

V 7.960±0.005 8.234±0.007 5 
AV 0.272±0.003 5 
B-V 0.815±0.003 1 
B-V 0.793±0.021 0.830±0.034 5 
V-I 0.83±0.02 1 
H p scatter 0.01 l a 1 
J 2M ass 6.585±0.024 6.746±0.030 6 
H 2M ass 6.218±0.046 6.362±0.034 6 
K 2M ass 6.151±0.026 6.267±0.017 6 
NUV magnitude 14.02 7 
FUV magnitude 7 

T eff (K) 5232 5073 5 

AT e ff(A — B) (K) 164±21 5 

log g 4.31 4.32 5 

lo S R HK -5.04 -5.07 3 

\ogR HK -5.07 8 

log^HK -5.06 -5.04 9 

usini (km/s) 0.6 1.3 3 

[Fe/H] 0.01 0.00 5 

A[Fe/H](^-B) 0.009±0.017 5 

Mass(M ) 0.91±0.03 b 0.88±0.03 b 4 

Age (Gyr) ~ 4 - 8 4 



References: 1 H ipp arcos ilPerrvman fc ESA] [1997); 2 
Ivan Leeuwenl (l2007Tl; 3 IDesidera et al.) d2006b[): 4 This Pap er; 
5 IDesidera et al.1 (l2006al); 6 2 MASS (iskrutskie eTall [2009); 7 
Galex dMartin et al.ll2005h : 8 iGrav et all lj2003T ): 9 lArriaeadal 
(I20T1 

a A+B 

b Deriv ed using the PARAM web interface (|da Silva et al.l 
[20061 ) 



orbit have an rms of about 6 m/s, and the periodogram 
does not show indication of additional significant periodici- 
ties (Fig. [2]). We then conclude that HD 106515A is orbited 
by a companion whose minimum mass is in the planetary 
regime. This should then be identified as HD 106515Ab. In 
the following section we will exploit the available astromet- 
ric data to further constrain the mass of the companion. 

The RV curve of the wide companion HD 106515B 
shows instead a scatter only slightly larger than internal 
errors (rms 8.7 m/s), without significant periodicities or 
long term trends (Fig. |3J). 



planet around HD106515A 
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Fig. 1. Upper panel: radial velocities of HD 106515A 
Overplotted the Keplerian best fit. Lower panel: Lomb- 
Scargle periodogram of the radial velocities. 

Table 2. Orbital parameters and results of fitting for RV 
of HD106515A. 



Parameter 


Our fit 


Mayor et al. 2011 


Period (d) 


3567±14 


3630 


K (m/s) 


160±3 


174 


e 


0.57±0.01 


0.60 


lo (deg) 


124±14 




TO (JD-2450000) 


1844±27 




msini (Mj) 


9.33±0.16 


10.50 


a (AU) 


4.43±0.01 




rms res (m/s) 


6.0 





5. Search for astrometric signature 

Considering the parameters from the RV orbit and the dis- 
tance to the system, the amplitude of the astrometric signa- 
ture for the RV minimum mass is about 1.2 mas. While we 
do not expect to be able to detect such astrometric motion 
in current data, it is possible to have a significant detection 
in case the orbit is seen close to pole-on and the actual mass 
is significantly larger than the projected mass. 

We considered all the relative astrometry measurements 
available in Washington Double Star catalog (kindly pro- 
vided by B. Mason), that span about 180 yr. Long period 
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Fig. 2. Upper panel: Residuals from best fit orbit vs time 
Lower panel: Lomb-Scargle periodogram of residuals from 
the Keplerian best fit of the radial velocities of HD1065f5A, 
with false alarm probability levels from bootstrap simula- 
tion overplotted. 



trends are clearly seen in both position angle and projected 
separation, as expected from the orbital motion of the wide 
pair (Fig. S]). The observed slopes will be used in Sect. [5] to 
constrain the binary orbit. Short term slopes as measured 
with Hipparcos are consistent with those based on the full 
dataset, supporting the lack of additional high amplitude 
astrometric perturbations on the timescales comparable to 
the mission lifetime. 

Residuals from long term quadratic slope show a much 
larger scatter before 1930 (Fig. [5]). We therefore considered 
only data taken after 1930 for the search of astrometric 
signatures due to the RV companion. Lomb Scargle peri- 
odograms of residuals in X and Y coordinates have very low 
power at period close to that of RV orbit. A possible period- 
icity is revealed at about 70 yr, especially in Y coordinate 
(the modulation can also be seen by eye in the projected 
separation vs time plot). Given the dense sampling of the 
astrometric measurements, this periodicity should not be 
related to the RV companion. 

If we limit our analysis to data from US NO (iFranzJ 

ll963t|Kallarakal et alJ[l969t Uosties et alJll978h and lHeintzl 
(|1987 ) (Fig. [7]) , to have a more homogeneous dataset and 
reduce the impact of absolute calibration errors that might 



Fig. 3. Upper panel: RV time serie of HD106515B Lower 
panel: Lomb-Scargle periodogram of RVs and false alarm 
probability levels from bootstrap simulation. 



be present when using different instrumentation, we have 
that the residuals have a dispersion of 13 mas in X and 20 
mas in Y, and of 14 and 15 mas in X and y respectively 
when taking the 70yr periodicity mentioned above into ac- 
count. These data cover 22 years. We then estimate an up- 
per limit on the astrometric amplitude with 10 yr period of 
about 30 mas, comparing the expected astrometric motion 
of the c ompa nion for different masses with the USNO and 
iHeintd (|1987l) dataset. This corresponds to a limit in mass 
of about 0.25 M . 

The tentative 70 yr periodicity with its possible am- 
plitude of about 50 mas would correspond to a very low 
mass star of about 0.1 M . The corresponding semimajor 
axis would be of about 17 AU, i.e. 0.5 arcsec on the sky. 
Our images taken with AdOptQTNG do not reveal stel- 
lar companions with masses larger than 0.15 — 0.2 M© at 
such a projected separation either around HD106515A or 
around HD106551B (Fig. [5$, therefore, the non-detection 
is not conclusive. The RV semiamplitude of a 0.1 M Q com- 
panion in 70yr orbit in edge-on circular orbit is about 700 
m/s, with RV difference that might reach 600 m/s over 10 
years. The real slope might be significantly smaller than 
this value, depending on the actual phase of the orbit at 
the time of the observations, inclination, eccentricity. The 
data of HD106515B does not reveal significant long term 
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Fig. 4. Upper panel: projected separation vs time for 
HD106515. Bottom panel: position angle vs time. In both 
panels filled circles are the measurements kept in the fit 
and the empty circles are those removed as being outliers 
in cither projected separation or position angle. In both 
panels, red continuous line is the linear fit imposing pas- 
sage through Hipparcos measurement, blue dotted line is 
the quadratic fit, purple dashed line is the slope as mea- 
sured from Hipparcos (baseline 3.25 yr; epoch 1991.25) ex- 
trapolated to the whole baseline of available data. 



trends (Fig. [3]), making very unlikely that such an addi- 
tional companion is orbiting around this star. The RV data 
for HD106515A are less adequate for the study of additional 
long term trends as the temporal baseline slightly exceeds 
one orbital period of HD106515Ab and the sampling after 
2009 is poor. Nevertheless, we do have specific indications 
of the presence of an additional long term trend in the data. 

We also note that no measurable differences are found 
in the stellar proper motion between the new reduc- 
tion of t he short-term Hipp arcos astrometric data for HD 
106515 ( van Leeuwen l2007t ) and the long-term Tycho-2 
data (|Hog et al.l l2000t ). The proper motion values in the 
two catalogs agree well with each other within the quoted 
uncertainties, thus no useful constraint can be obtained on 
the orbit and mass ratio based on the A/i technique (e.g. 
iMakarov fc Kar)lanll2005h . 



Fig. 5. Residuals from quadratic fit in Fig.0] Upper panel: 
projected separation, lower panel: position angle. 



6. Binary orbit 

To constrain the orbit of the HD106515 system, we con- 
sidered the long term relative astrometry (Sect. [3]) and the 
RV difference between the components. To derive this latter 
quantity, we measured the radial velocities of both compo- 
nents using the stellar template of HD106515A. Taking the 
orbital motion of the massive planet around the primary 
into account, such a difference results of ARVa-b = 739 
m/s. Internal errors derived from the error of the mean of 
the RV of the individual components and number of spec- 
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Fig. 6. Lomb-Scargle periodogram of the residuals from 
full-data quadratic fit in X and Y coordinates. Only data 
after 1930 were considered because of their higher precision. 
The short vertical lines at lOyr mark the RV period. 



tra and including uncertainty in the planet orbit are within 
10 m/s. Syst ematic effects due to spectral mismatch were 
estimated bv lNidever et al.l ()2002f ) to be of the order of 100 
m/s when using solar spectrum as template for the deriva- 
tion of the absolute RVs of FGK stars. This effect would be 
significantly smaller in our case thanks to the small tem- 
perature difference between the components and the very 
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Fig. 7. Residuals from full-data quadratic fit in X and Y 
coordinates. Only data from USNO (yearly averages) and 
Heintz 1987 are shown. 



similar metallicities and projected rotational velocities. The 
difference in the convetive blueshift and gravitational red- 
shift of the components am ount to just 24 m/s following 
Eq. 3 of iNidever et~aTI (|2002f ). Overall, the true error of our 
determination is likely within 50 m/s. 

Having the relative position and velocities on the plane 
of the sky and the relative velocity along the line of sight, 
we then derived the family of possible bound solutions as 
a function of the separation between the components along 
the lin e of sig ht (z) using the approa ch bv lHauser fc Marcvl 
(11999ft as in iDesidera etall (|201lh . The results are dis- 
played in Fig. IH1 

Critical semimajor axis for dynamical stability of 
planets around each co mponent, calculated following 
iHolman fc Wie^eTtl (fl999h . are shown in Fig. [TU] for the 
family of possible orbits shown in Fig. [SJ This quantity is 
larger than 40 AU for all the orbit£|. Therefore, the planet 
candidate is well within the stability boundaries for all pos- 
sible orbits of the wide binary. The potential astromctric 
candidate discussed in Sect. [5] is also well within the stabil- 
ity zone. However, if real, its presence should alter the RV 
difference, the astrometric trends and the mass of one com- 
ponent, with some effects on the family of possible binary 
orbits. 



7. Discussion 

The planet candidate around HD106515A, with its mini- 
mum mass of about 9.5 Mj, is one of the few with pro- 
jected masses close to deuterium burning limit. From avail- 
able data, we did not detect any additional planets in the 
system. The high eccentricity of its orbit and the solar-like 



1 At values of \z\ lar ger than 300 the large ecce ntricity of the 
binary orbit makes the|Holman fc WiegertT (Il999f ) equations no 
longer valid. 




1.0 1.5 

Projected separation (arcsee) 



Projected separation (AU) 

Fig. 8. Detection limits for stellar companions around 
the components of HD106515 from AdOpt@TNG images. 
Upper panel: AK vs projected separation in arcsec; Lower 
panel: Companion mass vs projected separation in AU. In 
both plots continuous line represent detection limits for 
HD106515A and dotted lines those for HD106515B 
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Fig. 9. Possible orbital parameters of the HD106515 system 
for various separations along the line of sight at present 
epoch. Top-left: semimajox axis; top-right: periastron of the 
orbit; bottom-left: eccentricity; bottom-right: inclination. 
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Fig. 10. Critical semimajor axis for dynamical stability for 
planets around the components of HD106515, for various 
separations along the line of sight at present epoch and 
corresponding binary orbit parameters. 



metallicity of its parent star are in agreement with the dif- 
ferences in the statisti cal properties of planets below and 
above 4 Mj found bv lRibas fc Miralda-Escudel (|2007Il 

The planet host has a stellar companion of similar mass, 
then HD106515Ab a dds to the growing census of exoplanets 
in mu ltiple systems (jDesidera k, Bar bicri 20(53; iRoell et ail 
I2012T) . The binary separation is quite wide and plausible or- 
bits leave dynamically stable zones up to 40-80 AU around 
the stars. This suggests a limited impact of the compan- 
ion on the planet properties but the moderate eccentricity 
might also be linked to Kozai interactions, which are effec- 
tive even for widely separated companion considering the 
old age of the system ([Takeda fc Rasiol2005l) . We also found 
from the analysis of the relative astrometry tentative indi- 
cation of an additional object with a period of 70 yr. From 
our data, we can not confirm the reality of this object, that 
might be a very low mass star, and infer the component 
around which it should be orbiting. 

The presence of a well separated companion with sim- 
ilar properties allowed us to perform a sensitive differ - 
ential abundance analysis (jDesidera et al.l 12004 l2006ajl . 
The lack of significant metallicity differences between the 
components extends the previous finding that large alter- 
ations of chemical abundances somewhat lin ked to the pres- 
ence of plan ets are not a common event (|Desidera et al.l 
2006a[ 1201 if) . There are 8 binary systems with giant plan- 
ets suitable for the comparison of che mical abundances (HD 
106515 and the 7 listed in Table 8 of lDesidera et alJl201lD . 
All of them have A [Fe/H] < 0.05 dex, which is significantly 
smaller than the typical difference in metallicity between gi- 
ant plan et hosts and nearby field s t ars (A [Fe/H] ~ 0.25 dex , 
see e.g. iFischer fe Valentil 120051: iNordstrom et al.l I2004D . 
This supports the primordial origin for the metallicity en- 
hanc ement of stars with giant planets (Fi scher &: Valentil 
1200511 . 

The binarity of HD106515 represents also a suitable op- 
portunity for the true mass determination HD106515Ab, 
thanks to the reference provided by HD106515B. The ex- 
pected astromctric amplitude is of about 1.2 mas for the 
minimum mass and about 10 mas at stellar/substellar 
boundary. 



From available relative astrometry, the orbital motion of 
the wide pair is clearly detected. From the analysis of resid- 
uals from the long term trend at epochs 1959-1980 where 
several high-quality data are available we put an upper limit 
to the mass of the companion of about 0.25 Mq. Much bet- 
ter astrometric precision can be obtained by more recent 
instrumentation. HD106515 is an i deal target for differen- 
tial astrometry u sing AO systems ( He lminiak et al.l [2009; 
IRoell et al.1l2010D and new interferometric instruments like 
PRIMA (|Quirrenbach et al.ll20Tll ). Subsequently, the com- 
bination of ground-based radial velocities and Gaia high- 
preci sion space-born e astrometric data might prove decisive 
(e.g.. ISozzettill20TTl and references therein). 

Therefore, there are very promising perspectives for a 
true mass determination of the companion of HD106515A 
in the coming years, then removing the ambiguity due to 
projection effects. The availability of true masses is rele- 
vant for a better understanding of the high-mass tail of the 
planetary mass function and the transition between planets 
and brown dwarfs. The direct detection of the companion 
is instead more challenging ev en for next generat ion planet 
finders as SPHERE or GPI (|Beuzit et alj 12010(1 . because 
of the small projected separations (< 0.2 arcsec) and faint 
luminosities implied by the old age of the system, unless 
the orbit is seen nearly pole-on and its mass is significantly 
larger than the minimum mass. 
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Table 3. Differential radial velocities of HD 106515 A 



HJD -2450000 


RV 


error 




m/s 


m/s 


1986.5429 


-198.8 


2.9 


2070.3741 


-180.9 


1.7 


2297.7099 


-135.2 


3.6 


2327.7449 


-116.3 


3.2 


2394.3984 


-94.2 


2.1 


2445.4047 


-91.1 


2.0 


2655.6473 


-40.4 


3.9 


2685.7174 


-42.8 


3.0 


2780.4680 


-30.7 


2.3 


3130.4615 


11.5 


2.9 


3370.7863 


24.5 


6.7 


3373.7381 


28.3 


2.5 


3454.5974 


39.8 


2.9 


3813.5819 


75.1 


4.7 


3898.4270 


75.5 


2.5 


4099.8020 


92.0 


2.1 


4160.6046 


95.3 


2.8 


4221.5300 


94.6 


2.1 


4251.4168 


97.6 


3.6 


4488.7506 


107.3 


3.2 


4545.4910 


105.4 


2.6 


4607.4763 


114.7 


2.8 


4609.4644 


119.3 


2.3 


4849.6782 


130.1 


3.0 


4962.4236 


121.2 


2.5 


5883.7810 


-109.7 


3.4 


5952.7662 


-92.9 


2.6 


5967.6966 


-98.9 


4.1 


6026.5722 


-100.2 


2.9 



Table 4. Differential radial velocities of HD 106515 B 



HJD -2450000 RV error 
m/s m/s 



1986.5544 


-19.0 


3.2 


2070.3879 


0.4 


1.9 


2297.7217 


2.6 


3.1 


2327.7578 


5.0 


3.5 


2394.5070 


5.5 


2.1 


2423.3945 


-0.6 


2.3 


2445.4169 


9.1 


2.8 


2655.6592 


0.6 


4.5 


2685.7307 


2.2 


3.3 


2780.4806 


-6.6 


2.6 


3130.4736 


6.8 


2.7 


3370.7985 


11.0 


4.9 


3373.7496 


-12.7 


3.8 


3454.6096 


-1.7 


3.7 


3813.5692 


1.5 


5.3 


3898.4402 


-0.1 


2.3 


4099.8161 


4.5 


2.3 


4160.6168 


-5.9 


2.9 


4221.5415 


-3.8 


2.4 


4251.4282 


-5.1 


4.5 


4488.7628 


7.6 


3.8 


4545.5030 


-12.1 


2.8 


4607.4885 


-9.3 


2.8 


4609.4782 


-4.3 


2.3 


4849.6896 


4.5 


3.1 


5952.7777 


17.2 


3.5 


5967.7085 


16.0 


5.1 


6026.5845 


-13.3 


3.4 



